Glueball and hybrid mass and decay with string tension below Casimir scahng 
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Lattice computations with excited SU(3) representations suggest that the confining gluon-gluon 
interaction complies with the Casimir scaling. The constituent gluon models have also been assuming 
the Casimir scaling. Nevertheless, inspired in type-II superconductors, we explore a new scenario 
for the gluon-gluon interaction where the adjoint string is replaced by a pair of fundamental strings, 
resulting in a factor of 2, smaller than 9/4. To test our proposal we construct a simple constituent 
gluon model, extrapolated from the funnel potential for quarkonium, and apply it to compute the 
wave-function of glueballs and of hybrid gluelumps. From the decay widths of quarkonium, we 
also extrapolate the decay widths of the glueballs. Our predictions apply to charmonia, lightonia, 
glueballs and hybrid gluelumps with large angular momentum. 



I. INTRODUCTION 



The glueballs are expected from QCD, and indeed 
they are observed in pure gauge lattice QCD simula- 
tions. The first experimental evidence for glueballs is 
an indirect one, the pomeron which explains the high 
energy scattering of hadrons, 0, H H g. Systematic 
comparisons of the glueball masses computed in Lattice 
QCD with the pomeron trajectory, initiated by Llanes- 
Estrada, Cotanch, Bicudo, Ribeiro and Szczepaniak, con- 
firm that the lattice QCD glueballs comply with the 
pomeron 0, 0, 0, ■ Direct experimental evidence for 
glueballs is still controversial , nevertheless in few years 
two major collaborations will update the experimental 
search for glueballs, PANDA at GSI and GLUEX at 
JLAB. Different theoretical approaches to glueballs are 
trying to match the experimental effort. 

Here we aim at a simple and plausible model for glue- 
balls and for other hadrons, with large angular excita- 
tions. Many of these hadrons remain to be observed, 
see Table 13 We first review the possible glueball mod- 
els, extending the string mechanism for quark-antiquark 
confinement. The quark model describes phenomenol- 
ogy, simply assuming constituent quarks, a short range 
coulomb interaction and a long range confining string in- 
teraction. This model is depicted in Fig. ^a). Nambu 
and Jona-Lasinio showed that the massless fermions 
can have a mass generated by the spontaneous breaking 
of chiral symmetry, with the Schwinger-Dyson technique, 
equivalent to a second-quantized canonical transforma- 
tion Indeed this has been successfully appHed 
to quark models 0, 0, l^l , suggesting that quark mod- 
els may comply both with the PCAC theorems and with 
the success of quark model spectroscopy 0, Q, 0, . 
Moreover the quark-antiquark confining string, consist- 
ing of a color-electric flux tube was suggested by Nielsen 
and Olesen ^^l- Nielsen and Olesen conjectured that 
the magnetic flux tube vortex of type-II superconductors 
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FIG. 1: Basic models b) c) and d) for glueballs, extended 
from the string model a) for quark-antiquark confinement in 
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could be extended to color-electric flux tube strings in 
QCD. 

Inspired in the constituent string-confined quark model 
a), three models for glueballs are plausible. They are re- 
spectively depicted in Fig. ^ In model 6), a constituent 
gluon pair is confined with an adjoint string, in model c), 
a constituent gluon pair is confined with a pair of fun- 
damental strings and in model d), a closed fundamental 
string loop is rotating , without any constituent gluon 
at all. Models b) 0, |i. El El HO] , with a pair of con- 
stituent gluons confined by an open adjoint string, and 
d) 0, 123I , with a closed fundamental string and no con- 
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TABLE I: Experimental (Review of Particle Physics) avail- 
able masses M and widths F, for charmonia, lightonia, K* , 
glueballs and gluelumps, with maximal s and j = I -\- s. 



states / J^*^ 


0++ 1— 


2++ 


3 — 


4++ 5— 


(^Vi annrmiiiTTi A/T 


3096 


3556 






Charmonium, P 


0.1 


2.1 






lightonium 1=0, M 


783 


1275 


1667 


2034 


lightonium 1=0, P 


8.5 


185 


168 


222 


lightonium 1=1, M 


775 


1318 


1689 


2010 


lightonium 1=1, P 


150 


107 


161 


353 


lightonium K* , M 


892 


1430 


1776 


2045 


lightonium K* , P 


51 


105 


159 


198 


Glueball, M 










Glueball, P 










Gluelump, M 










Gluelump, P 











stituent gluon, and also hybrid models HSI^Ii have been 
studied in the literature. The model c) is explored in this 
paper. 

An important ingredient of models h) and c) is the mas- 
sive constituent gluon. There is evidence for a massive- 
like dispersion relation for the gluon, with a mass rang- 
ing from 700 MeV to 1000 MeV, both from laticce QCD 
[25ll2^|2? . 28] and from Schwinger-Dyson equations |29l| . 
This suggests that the Anderson-Higgs mechanism is pro- 
viding a consistent framework for the mass generation of 
gauge bosons. For instance in the Meissner effect of su- 
perconductors, the photon has a mass. It is plausible that 
simple non-relativistic constituent gluon models may in- 
deed describe glueballs. 

Another ingredient, necessary for model h) only, is the 
adjoint string. The adjoint string is a natural extension 
of the fundamental string assumed in models a), c), d). 
It consists of an excited string with a colour-electric flux 
larger than the one of the fundamental string. Two pos- 
sible contributions for colour-colour interaction are the 
shorter range one gluon exchange interaction, and the 
longer range confining string interaction, proportional to 
the colour-electromagnetic energy density. Both are pro- 
portional to A • A, one of the two QCD Casimir invari- 
ants. Moreover Bali 's^, |^ showed, in pure-gauge lat- 
tice QCD, extending the lattice gauge field links from the 
fundamental representation to several representations of 
SU(3), that the colour-colour interactions comply with 
the Casimir scaling, they are all proportional to A • A, 
except for smaller than 5% numerical errors. 

Nevertheless we propose here a new model c) with the 
same pair of constitue nt g luons of b) , but with a gluon- 
gluon potential lower |32l | than the Casimir scaling of 
model c). Importantly, model c) is not excluded by the 
lattice QCD computations published so far. The Casimir 
scaling for the octet-octet interaction was observed in lat- 
tice discretizations of QCD where the fundamental links 



are replaced by adjoint links, 8x8 matrices of the adjoint 
representation of SU{i). Therefore these lattice compu- 
tations are not adequate to simulate our model c) which 
assumes fundamental strings only. 

Like the seminal paper of Nielsen and Olesen [Trj . 
model c) is inspired in type-II superconductors |3^ . In 
type-II superconductors the magnetic fiux is allowed to 
penetrate, and it is localized in vortices, with flux 

4>n = , (1) 

where the number n quantifies the magnetic flux. Rather 
than having a single vortex with a high flux, it is ener- 
getically favourable to have several elementary vortices, 
say a square lattice of vortices, each with the minimal 
flux (f>i . This occurs because the energy of each vortex 
is proportional to the square of the magnetic flux. For 
instance a vortex with n = 2 does cost the double of the 
energy of two elementary vortices with n — I, while both 
scenarios have the same flux (/)2 = 2^i. Back to our con- 
stituent gluon potential, it is energetically favourable to 
have two fundamental strings, with scaling factor of 2, 
than one adjoint (excited) string, with a scaling factor of 
9/4 = 2.25. 

In the remaining of this paper we specialize in the pro- 
posed model c), extending the basic constituent quark 
model, from the quark-antiquark meson, to the gluon- 
gluon glueball and to the quark-antiquark-gluon hybrid. 
Because we also study light mesons, our kinetic energy 
is relativistic. We specialize in systems with maximal 
angular momentum, because the string length is large 
and the details of the conflning string become relevant. 
These large J states also allows us to avoid the details 
of spin tensor potentials. The experimentally confirmed 
resonances only go up to J = 4 for light mesons and 
J = 2 for charmonium, see Table HJ therefore most of our 
results are predictions. 

In Section II we discuss the parameters of our model. 
The scaling of 2, the experimental meson masses, the 
glueball masses an parameters obtained in the lattice 
and in the Schwinger-Dyson equations, together with the 
pomeron intercept determine all the parameters of the 
model. In Section III we compute the masses (up to 5 
GeV, close to the energy limit of the future experiments) 
and the mean radius of charmonia, lightonia, glueballs 
and hybrids. We also estimate decay constants in Sec- 
tion IV. In Section V we conclude. 



II. THE PARAMETERS OF THIS 
CONSTITUENT QUARK AND GLUON MODEL 

A. Mesons 

Our starting point to parametrize the potentials and 
masses is the funnel potential for charmonium 34j. The 
charmonium spectrum and wavefunctions are usually de- 
termined solving the non-relativistic radial Schrodinger 
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equation for two-body systems, 



Mtc' 



+ 



h"^ 1(1 + 1) 



Vni{r) = rRni{r) , Mt = ml + m2 , fi 



+ V{r) Vnl{r) = EnlVnl 

nil "^2 



the Chew-Frautschi plot. A constant C term in the po- 
tential is also needed to fit the spectrum. Here the C 
sums the effect of spin-tensor potentials and of vacuum 



mi + 7712 



However we are also interested high angular excitations, 
and therefore it is more convenient to upgrade the kinetic 
energy to the relativistic one. In the center of mass frame 
this amounts to replace in eq. 



Mtc' 



-f 5 > ^mlc^ + p2 c2 4- ml + 



p2 



(3) 



Again we can perform the standard angular separation 
of the Schrodinger equation, where the squared momen- 
tum is separated in a simple second derivative and in a 
centrifugal barrier. 



r r 



lil + l) 



(4) 

thus the equation remains a radial equation, the spher- 
ical harmonics can be factorized. The arguments of the 
square roots in eq. © are positive definite and therefore 
there is no technical difficulty in computing the relativis- 
tic kinetic energy (|3Jl. A possible method of finding a 
function f{M) of a matrix M uses the eigenvalues Xi 
and eigenvectors. 



/(Af) = ^|z)/(A,)(* 



(5) 



Thus, once the squared momentum in eq. Q is diagonal- 
ized, either with spherical Bessel functions or with finite 
differences, the relativistic kinetic energy operator of eq. 
0, a function of the squared momentum, can be com- 
puted with the method of eq. . Then we finally solve 
the schrodinger equation, diagonalizing the hamiltonian. 
Numerically the relativistic code is only twice as slow as 
the non-relativistic one, because we perform two diago- 
nalizations. Importantly, the relativistic hamiltonian is 
adequate to stydy lighter systems, with light quarks or 
gluons. 

The parameters of this model include the charm mi 
and anti-charm m2 identical masses, and a Coulomb po- 
tential, a linear potential, and a constant potential, 



Vir) 



hcrr + C 

r 



(6) 



The Coulomb potential — — can be extracted from the 
perturbative QCD one gluon exchange potential and 
from the non-pertubative Luscher fluctuations of the 
string [35|. All Coulomb contributions are summed in 
the parameter a. The linear potential is parametrized 
by the string tension a. Both can be extracted from 
lattice QCD. Moreover the linear potential needs some 
Coulomb interaction to linearize the Regge trajectory in 



effects. Moreover, in Lattice QCD the constant term in 
,tJae potential is not determined, and in a sense this con- 
^stitutes a free parameter fo the funnel potential. Here 
we use the parameters mi — m2, a, a and obtained by 
Lucha and Schoberl 34] to fit the the charmonium spec- 
trum. Because our kinetic energy is relativistic, we refit 
the constant C . Our charmonium parameters are show 
in Table in. 

Our next class of hadrons is the lightonia, mesons com- 
posed of light quarks. In this case the masses are ex- 
pected to be slightly larger than the third of the lighter 
baryon masses, see for instance the modelling of baryons 
by Isgur and Karl j36l |. In what concerns the Coulomb 
and string tension they may be slightly different than the 
ones of Charmonium. The parameters that we used to 
fit the experimental spectrum are again show in Table 
im While the same Charmonium string tension can be 
used, a larger Coulomb potential is needed to linearize 
the light meson Regge trajectory, and the constant po- 
tential is also slightly refitted. This is consistent with 
the Fermi-Breit short range Coulomb potential, which 
increases for light masses. 



B. Glueballs and hybrids 

Our next class of hadrons is the one of glueballs. The 
mass of the constituent gluon is determined with the 
Schwinger-Dyson equation for the gluon propagator, [2fll | 
equivalent to the Bogoliubov-Valation mass gap equation 
|37| for the gluon propagator, and with lattice QCD com- 
putations of the gluon propagator [23, H^- Althought 
the gluon propagator is not gauge invariant, and it is not 
clear whether gauge symmetry is broken or not, there is 
evidence for a mass ranging from 700 MeV to 1 GeV in 
the dispersion relation of the gluon. This scale is also 
present in the lattice QCD determinations of the first ex- 
citation of the quark-antiquark or three-quark potentials 
[3M Isgf . Here we use a mass of 800 MeV, lighter than 
the mass of a light quark-antiquark pair. In this case the 
constituent gluon might be relatively stable against the 
decay to a quark-antiquark pair. 

To address the gluon-gluon potential, we first review 
the Casimir scaling. In the fundamental representation 
of SU(3), acting in the quark colour triplet, the Lie Al- 
gebra is represented by the 3x3 Gell-Mann fundamental 
matrices A". The product of the fundamental A" is. 



a \ h 



A" A 



(7) 



where the antisymmetric structure constants ^ ^ are an- 
tisymmetric and independent of the representation. The 
constant tensors 1(5" and d°-^^ are symmetric and de- 
pend on the representation. Another relevant represen- 
tation is the adjoint representation. In the adjoint repre- 
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sentation, acting in gluons, the Lie Algebra matrices are 
represented by the 8x8 adjoint matrices, 

With eq. iQ, and using the properties of the antisym- 
metric f'^'"^ and of the symmetric d° '"^ we find that the 
Casimir invariant in the fundamental and adjoint repre- 
sentation are respectively. 



^adjoint ^adjoint ^ ' V^/ 

and they differ by a factor of 9/4, the Casimir scal ing, . 
The Casimir scaling was indeed observed by Bali [sj 
in lattice QCD. Bali showed that the energy density is 
proportional to the A°A". This Casimir Scaling is the 
simplest possible result, and Semay 40] showed that it 
occurs if the color-electric flux-tube thickness is essen- 
tially the same for all possible quantized fluxes. 

Nevertheless two fundamental strings cost less energy 
than one adjoint string. The scaling of the present model 
is 2 < 9/4. Because the mass of our constituent gluon 
is between the mass of the charm and light quarks, and 
the charmonia and lightonia used the same string ten- 
sion, for the string tension of the glueballs we double the 
parameter a already used for the quarks. 

The parameter a of the glueball Coulomb potential 
can also be indirectly estimated from the masses of glue- 
balls in lattice QCD. The Lattice QCD result for the 
hyperfine splitting between the masses of the 2++ and 
O"*""*" is of the order of 0.7 GeV. This splitting can be 
compared to the hyperfine splitting between the 1 
and ^ mesons, which is of the order of 0.4 GcV for 
light mesons (once chiral effects are subtracted) and is 
of the order of 0.1 GeV for charmonium. Using the 
one gluon exchange hyperfine splitting potential propor- 
tional to fiAii^l^W 5^. . 5^, similar to the one derived for 
mesons |41|| . Cornwall and Soni and Kaidalov and Si- 
monov 0,1421, found an hyperfine splitting, between the 
2++ and 0++ glueballs, of the order of 1.0 Gev. Assuming 
that the hyperfine splitting potential is essentially due to 
the pertubative one gluon exchange potential, this sug- 
gests that the Coulomb potential should be suppressed 
by 30%. However, with the scaling factor of 2 , rather 
than the Casimir Scaling of 2.25, we already decrease the 
hyperfine potential by a factor of more than 10%. More- 
over these authors used gluon masses of the order of 0.5 
to 0.6 GeV. Our gluon masses are larger (because they 
are estimated in Lattice QCD and in Schwinger Dyson 
truncated QCD equations) and this further decreases the 
hyperfine potential. Thus we are confident that a glue- 
ball Coulomb potential with 2 x the same a already used 
in charmonium is consistent with a correct glueball hy- 
perfine splitting. 

In what concerns the constant potential C, it is fitted 
with the intercept of the glueball trajectory, which is ex- 
pected to reproduce the pomcron. The equation for the 



TABLE II: Parameters of the constituent model for glueballs, 
charmonium and gluelumps. 



parameters 


mi [GeV] 


ma [GeV] 


a 


a [GeV^] 


C [GeV] 


charmonium 


1.64 


1.64 


0.27 


0.25 


-0.69 


lightonium 


0.44 


0.44 


0.35 


0.25 


-0.80 


glueball 


0.80 


0.80 


2x0.27 


2x0.25 


-0.52 


gluelump 


0.80 


3.28 


2x0.27 


2x0.25 


-1.59 



J 




FIG. 2: Chew-Frautschi plots of the leading trajectories for 
charmonia (black), lightonia (dark), glueballs (grey) and hy- 
brid gluelumps (light). 

pomcron trajectory in the J, t — AP space is 0,0], 

J = ap{t), (10) 
ap(t) = 1.08 + 0. 25i . (11) 

The intercept ap (t) is of the order of 1 , and this explains 
the high energy hadronic cross sections. In particular we 
may ignore the small decimal digits .08 because they may 
be due to double pomcron exchange. The pomcron is also 
expected to correspond to a series of glueball masses. 
Therefore our final parameter C is fitted to the intercept 
ap{0) = 1. Again our glueball parameters are shown in 

Table nn 

In what concerns the hybrids, it is possible that we 
have a gluelump. The charm and anticharm pair may 
form a colour octet with a mass of 3.3 GeV, quite 
larger that the constituent gluon mass. Then the charm- 
anticharm octet is essentially stopped at the centre of the 
mass of the hybrid system, with the glueball orbiting it, 
attached by a double fundamental string. Such a system 
of a gluon attached to a heavy colour octet is called a 
gluelump Theoretically, this is a very interesting 

object because it is even simpler than a Glueball. All 
the parameters of the gluelump potential are the same of 
the glueball potential, except for the constant potential 
C witch may be refitted. 

The constant potential C is fitted to produce a split- 
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ting of 1 GeV between the first gluelump state and the 
first S — 1 charmonium state. This spUtting has been ob- 
served in lattice QCD for quenched quark potentials [38l |. 
Since in the gluelump we ignore the motion of the charm 
and the anticharm, the constant C not only accounts for 
the spin-tensor interaction, and for vacuum effects, it also 
accounts for the kinetic and potential energy internal to 
the cc pair. Again our gluelump parameters are shown 
in Table in 



III. MASSES AND RADII OF GLUEBALLS, 
HYBRIDS AND QUARKONIA WITH LARGE J 



in particular the slope is close to the slope of 0.25 of the 
soft pomeron model of Donnachie and Landshoff [43l| . 
The model b) with Casimir scaling would produce the 
trajectory, 

j ^ 1.00 + 0.20t . (13) 

For the same gluon masses of 0.8 GeV, model b) with 
Casimir scaling has a smaller slope, smaller than the Don- 
nachie slope. 

The four different trajectories, respectively for the 
charmonium, lightonium, glueball and hybrid are also 
displayed in Fig. |21 



We solve the radial Schrodinger equation, with rel- 
ativistic kinetic energies Q, with a finite difference 
method, which transforms the differential equation in 
a simple matrix eigenvalue Method. The parameters 
of our constituent model for glueballs, charmonium and 
gluelumps are show in Table Hll 

We are particularly interested in large angular mo- 
mentum states, with maximal S = Si + S2 and with 
J = L + S . The maximal spin is 5 = 1 for the lighto- 
nium and charmonium. The maximal spin is S* = 2 for 
the glueball and for the gluelump. We do not compute 
other spin-orbit combinations or the radial excitations, 
because they belong to daughter trajectories, with larger 
decay widths, which are harder to identify experimen- 
tally. 

The possible quantum numbers of the glueball are fur- 
ther constrained, because the gluons are bosons. Since 
in a colour singlet the colour wavefunction is symmetric, 
and because we are only considering here the maximal 
S = 2 which is also symmetric, the angular wavefunction 
needs also to be symmetric. Therefore the maximal J 
trajectory starts at J = 2(L = 0) and continues with 
J = 4(L = 2), J = 6(L = 4)--- 

In what concerns the gluelump, the spin of the charm 
and anticharm pair is not expected to affect the gluelump 
mass significantly. For a maximal S we consider that the 
cc pair has spin 1 and that the total spin of the ccg system 
is 2. Then the leading trajectory has J ~ 2{L = 0), J ~ 
3(L = 1), J = 4(i = 2) - •• 

We compute the energies Eqi (up to the first state 
above 5 GeV) and the corresponding wave-functions 
Roi{r). The energies and the radius mean square are 
displayed in Tables UTTI and Hvl All our units are in pow- 
ers of GeV, for h — c ~ 1. For instance r = lQGeV~^ = 
1.97/m. 

We also compare our glueball trajectory with the one 
with the Casimir scaling. The potential C is refitted in 
the Casimir scaling case, to maintain the intercept equal 
to 1. The new constant potential is C = —0.39, and only 
the slope of the Chew-Frautschi trajectory changes. Our 
model c), with a pair of fundamental strings, yields the 
trajectory 

J = 1.00-|-0.24i , (12) 



IV. ESTIMATING DECAY WIDTHS OF 
MESONS AND GLUEBALLS WITH LARGE J 

We now estimate the decay width of the predicted 
hadrons, and this includes several hadrons with large J 
which have not yet been detected. Let us consider for 
instance the decay of a glueball 0, 0|. We may as- 
sume that the decay is initiated when, either a gluon in 
the colour-electric flux tube, or a "massive" constituent 
gluon, is transformed into a quark-antiquark pair with 
the QCD quark-gluon coupling. 

Then the glueball is transformed into an excited meson 
or a hybrid. Because this is not a stable system, it also 
decays, again with the same quark-antiquark pair cre- 
ation. Essentially we have a cascade of decays, because 
the first decay products are unstable systems. This ac- 
counts for the final decay product of several pseudoscalar 
mesons 

Let us consider the decay scenario where string break- 
ing dominates over the direct constituent gluon decay. 
If the constituent gluon mass is lower than the con- 
stituent quark-antiquark mass, such a dominance is pos- 
sible. Here rrig = 0.8 GeV < 2 m, = 0.88 GeV. While 
the problem of the constituent gluon decay remains to 

Ro,6(r) [Gev3/2] 
.09 



0.06 
. 03 




FIG. 3: Example of the radial part Rni(r) of the wave- 
function for a glueball with a high angular excitation (1=6). 
The total wave-function is obtained multiplying the radial 
wave-function by the spherical harmonic Yl^{6, (f>) 
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be solved, here we neglect this effect. Therefore our re- 
sults for the decay widths are lower bounds on the actual 
decay width. 

In this scenario a string is cut and a qq pair is cre- 
ated. Each quark remains at one end of the open parts 
of the string. Statistically, it is plausible that the prob- 
ability to produce a qq pair is proportional to the string 
length. It occurs that 2x the string length in the glue- 
ball g g, is close to the string length in the q q case, with 
the same angular momentum. Thus, in the glueball with 
a double fundamental string, the total string length is 
similar to the string length of the light meson with the 
same angular momentum. So the probability to have 
the breaking of a string is similar in these glueballs with 
J^'-^ = 2++, 4++, 6++... is respectively similar to the 
probability of having a string breaking in light mesons 
with JP'^ = 1 — , 3—, 5~".... 

To compute the decay widths, one would still have to 
evaluate the overlaps with the decay products, and the 
available phase space. In multiple particle products, it 
remains difficult, in the present state of the art of glue- 
balls, to study quantitatively the decays. Therefore we 
will essentially focus in the relation of the string length 
with the decay width. 

For an educated guess it is convenient to study the de- 
cays of conventional mesons. The decay widths extracted 
from the Review of Particle Physics 'iq are show in the 
TableQ] The corresponding string lenghts can be approxi- 
mated by the RMS of Table llVl To show the dependence 
of the decay width in the string length, a plot is drawn 
in Fig. 21 

In Fig. ^it is clear that the decay width of the light 
mesons is well fitted by a linear increase with the string 
length. Fitting the average of the different / = 0, / = 1 
and K* lightonium decay widths we get 

r = -f{RMS - ro) ± AF , (14) 



r [GeV] 



FIG. 4; Decay widtli of angular excited resonances as a func- 
tion of the theoretically estimated RMS. The decay widths 
are extracted from the Review of Particle Physics llq. The 
black dots are the average of the different J = 0, / = 1 and 
K* lightonium decay widths. The vertical lines span from the 
maximum to the minimum of these lightonium decay widths. 
The dashed line corresponds to the best linear fit of the light 
data. The gray dots represent the charmonium decay widths. 

lightonia and glueballs with large angular excitations are 
of the order of 10% of the respective masses. 

However, in what concerns charmonium, the D mesons 
are quite heavier than the light mesons. The phase space 
of the charmonium decay products is quite different from 
the phase space of the lightonium decay products, and 
we should not apply this simple linear rule, with the 
same parameters, to the charmonium decays and to the 
gluelump decays. Moreover, only the \~~ and 2++ char- 
monium decay widths are known, and this is not sufficient 
for a precise linear fit of the decay widths. 



where 7 = 0.05 ± 0.01 GeV^ and ro = 1.4 ± 0.6 GeV"!. 
The error in the parameters 7 and tq are statistical errors. 
The error AF = O.lGeT^ includes the details of the decay 
processes. 

If we extrapolate this linear growth to mesons with 
higher angular excitations, the decay width grows up to 
0.68 GeV for the J^'^=15— meson made of light quarks, 
with a RMS of 15.24 GeV"!. The decay widths of the 
light mesons can also be extrapolated to glueballs. In 
this case, because the glueball has a double fundamen- 
tal string, for a similar RMS, the glueball has a string 
lenght twice as long. Thus, for glueballs, the input vari- 
able RMS of eq. (|14|l must be doubled. For instance, 
decay widths respectively of 0.14, 0.33 and 0.48 GeV are 
expected for the first three glueballs of our series, provid- 
ing the string- breaking decay process is dominant. If the 
direct decay of constituent gluons turns out to be impor- 
tant, the glueball decay widths are larger. Our predicted 
widths for lightonia and for glueballs are shown in Table 
fVl Essentially, we predict that the decay widths of both 



V. CONCLUSION 

With a simple constituent model, we predict the 
masses and mean radius of charmonia, lightonia, glue- 
balls and gluelumps. We We also predict the decay 
widths of lightonia and glueballs. The ingredients are 
constituent quarks and gluons, and the fundamental 
SU(3) chromo-electric flux tube. 

We compute the spectrum and wave- functions of 
hadrons with angular excitations and maximal J. These 
hadrons with no radial excitations may be easier to iden- 
tify in the lattice and in the experiments. Because we 
construct all the hadrons with the same confining strings, 
we are able to extrapolate the decay widths of mesons to 
glueballs, in a unified framework. This model can also 
be applied to three-gluon glueballs ^-^id other many- 
particle systems. 

Essentially we expect the angularly excited lightonia 
and glueballs to have decay widths of the order of 10% 
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of the respective masses. The largest decay width we 
get, for hadronic masses up to 5 GeV, are of the order 
of r = 0.6 GeV. Notice that, because we neglected the 
direct decay of the constituent gluons, the actual decay 
widths of glueballs and hybrids may be larger than the 
ones predicted here. We nevertheless point out that the 
decay width F = 140 MeV of our 2++ glueball is identi- 
cal to the one predicted by Cotanch and Williams with 
vector meson dominance jSfl, |5l| • 

However the large mass of the glueballs and of the 
other studied hadrons enable final decay products into a 
large number of mesons. The first decay is expected to 
result in excited hadrons, and a cascade of decays into 
lower states is expected to produce several pions in the fi- 
nal state. Because the partial wave analysis may be quite 
difficult when, say, more than four pions are produced, 
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TABLE III: Predicted energies in GeV of charmonia, lightonia, glueballs and gluelumps as a function of J. We limit our study 
to maximal J states (no radial excitations and J = L + Si +52) and to energies up to 5 GeV. 



hadron \ J 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


14 


15 


Charmonium 


3.12 


3.53 


3.87 


4.1 


4.46 


4.74 


4.88 


5.10 
















lightonium 


0.77 


1.27 


1.68 


2.05 


2.40 


2.72 


3.04 


3.34 


3.63 


3.91 


4.18 


4.45 


4.70 


4.96 


5.21 


Glueball 




2.08 




3.54 




4.67 




5.68 
















Gluelump 




4.12 


4.88 


5.47 

























TABLE IV: Predicted RMS — y (r^) in GeV ^ ~ 0.20 fm of charmonia, lightonia, glueballs and gluelumps as a function of 
J. 



hadron \ J 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


14 


15 


Charmonium 


2.24 


3.35 


4.34 


5.26 


6.13 


6.96 


7.77 


8.54 
















lightonium 


2.73 


4.06 


5.22 


6.28 


7.27 


8.21 


9.11 


9.96 


10.79 


11.59 


12.37 


13.11 


13.83 


14.55 


15.24 


Glueball 




2.10 




4.06 




5.68 




7.13 
















Gluelump 




1.86 


2.82 


3.64 

























TABLE V: Estimated decay widths in GeV of lightonia and glueballs as a function of J. Only the string breaking effect is 
considered. The statistical and systematic error bars are detailed in eq. H14|l . 



hadron \ J 





1 2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


14 


15 


lightonium 




0.06 0.13 


0.19 


0.24 


0.29 


0.34 


0.38 


0.42 


0.46 


0.50 


0.54 


0.58 


0.61 


0.65 


0.68 


Glueball 




0.14 




0.33 




0.49 




0.64 

















